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We report on reproducible shock acceleration from irradiation of a A = 10 /xm CO 2 laser on 
optically shaped H 2 and He gas targets. A low energy laser prepulse (/ < 10^^Wcm“^) was 
used to drive a blast wave inside the gas target, creating a steepened, variable density gradient. 

This was followed, after 25 ns, by a high intensity laser pulse (/ > lO^^Wcm”^) that produces 
an electrostatic collisionless shock. Upstream ions were accelerated for a narrow range of prepulse 
energies (> 110 mJ & < 220mJ). For long density gradients (> 40/xm), broadband beams of He^ and 
were routinely produced, whilst for shorter gradients (< 20/im), quasimonoenergetic acceleration 
of proton was observed. These measurements indicate that the properties of the accelerating shock 
and the resultant ion energy distribution, in particular the production of narrow energy spread 
beams, is highly dependent on the plasma density profile. These findings are corroborated by 2D 
PIC simulations. 


The light and thermal pressure associated with intense 
lasers causes compression and heating when incident on 
an overdense plasma. This can generate a piston that 
launches electrostatic collisionless shocks into the plasma 
[1]. Upstream ions can be reflected off the shock creating 
an ion population accelerated to twice the velocity of the 
driving shock PHS]. This mechanism has been demon¬ 
strated in recent experiments using intense CO 2 lasers 
interacting with gas jets laiTi where multi-MeV proton 
beams with energy spread smaller than 4% have been 
reported. 

Gas jets have also been proposed as a source of high 
purity, high-Z ion beams. This is in contrast with the 
multiple species beams generated from solid targets via 
sheath acceleration m- For solid targets, protons from 
surface contaminants are preferentially accelerated due to 
their higher charge-to-mass ratio, making the production 
of impurity free high-Z beams challenging. Generating 
high purity helium ion beams from gas jets would be of 
interest for nuclear and medical physics applications [TO] . 

Laser driven longitudinal m and transverse [siiiaiis] 
acceleration of helium ions from laser irradiation of gas 
jets has been observed experimentally. However, these 
previous studies were conducted in underdense plasmas 
requiring laser intensities > lO^^Wcm”^. Gollisionless 
shock acceleration, on the other hand, allows the produc¬ 
tion of directed beams at lower intensities with sizeable 
number [6], and with favourable intensity scaling [5]. 

It is thought that the generation and subsequent prop¬ 
erties of shock accelerated ion beams are highly depen¬ 
dent on the initial plasma density distribution BM- 
Gas jet targets used for laser-plasma interactions typi¬ 
cally have initial density scale lengths / > 100 /xm, which 
is unsuitable for efficient shock generation. In previous 


experiments using intense GO 2 lasers, this profile was 
modified by the pulse train inherent to the laser system 
BB, which is difficult to control and reproduce. Re¬ 
cent breakthroughs in GO 2 laser technology have made 
single pulses of intense radiation m possible. This al¬ 
lows reproducible interaction of the laser with the target. 
However it requires a different method to modify the den¬ 
sity profiles. One method investigated previously is the 
use of optically generated hydrodynamic blast waves via 
laser solid interaction [liiiT]. 

In this paper, we demonstrate the use of a controlled 
low energy prepulse focused in dense gas to generate blast 
waves that hydrodynamically shape the gas targets to 
make them suitable for ion acceleration. Interaction of 
the high intensity laser pulse with the modified plasma 
density profile leads to reproducible shock acceleration. 
In particular, we present the first measurements of shock 
accelerated, > 1 MeV, helium ions. Using hydrogen, we 
also demonstrate that shock acceleration can lead to both 
broadband and quasi-monoenergetic ions beams depend¬ 
ing on the steepness of the initial plasma density gradient 
The role of the density gradients is further investigated 
with PIG simulations. 

The experiment was performed at the Accelerator Test 
Facility at the Brookhaven National Laboratory. A lin¬ 
early polarized GO 2 laser beam, A = 10.3 /xm, with 
energy up to IIJ per pulses of 5ps full-width-half¬ 
maximum (FWHM), provided a peak power of 2.2TW. 
An //3 parabolic mirror focussed the laser to a spot wq = 
65 /xm, resulting in peak intensity I = 2.5 x 10^^ Wcm“^ 
(ao = 1.4). The laser was focussed at the center of a 
helium or hydrogen gas-jet from a 1 mm diameter super¬ 
sonic nozzle, ^ 700 /xm above the the nozzle. At this po¬ 
sition the initial gas density profile, along the laser axis. 
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was near triangular with a linear gradient of ^ 1 mm up 
to a peak fully-ionized plasma density of rzg ^ 2 nc m- 

The ion beam was characterized using a laser-axis 
Thomson parabola spectrometer, opening angle of 2 x 
10 “^ sr, and detected with BC-408 polyvinyl-toluene 
scintillator coupled to an EMC CD camera. The ion-to- 
detected-photon yield was calibrated using the Tandem 
proton accelerator at Stony Brook University [19]. The 
calibration was scaled to He+ by introducing a scaling 
factor of 0.3 ±0.1 [20] . 

A laser prepulse collinear to the drive pulse was gen¬ 
erated to shape the gas. It arrived 25 ns before the 
main laser pulse and contained a variable energy, Epp^ 
up to IJ in 5ps. The resulting plasma density distri¬ 
bution was characterized using optical probing. Two, 
time-adjust able, ^ 10 ps long, A = 532 nm synchronized 
pulses were used for shadowgraphy and interferometry. 
They were set to allow characterization just before and 
after the drive pulse on a single shot. 

Fig# shows the experimentally observed plasma den¬ 
sity map 100 ps before the main laser plasma interaction 
(LPI) using He, similar images are found using H 2 . For 
this shot, a prepulse Epp ^ 220 mJ was incident on he¬ 
lium gas of peak neutral density uhc = 0.8 x 10 ^^ cm“^, 
corresponding to a fully ionized plasma density Ue = 
1.6 Uc- The plasma density was obtained from interfer¬ 
ometry assuming cylindrical symmetry around the laser 
axis, thus allowing Abel inversion. The optical path dif¬ 
ference due to the unionized gas is small compared to 
that of the plasma and is therefore neglected. Plasma 
formation was observed over a large volume with a peak 
density significantly lower than the fully ionized density, 
and shows no structure. However, probing the plasma 
300 ps after the LPI revealed the characteristic high den¬ 
sity shell of a blast wave (fig. &)• The shell is observed 
to expand at velocity < 10^ms“^. Hence the dramatic 
change in the plasma appearance cannot be due to hy¬ 
drodynamic motion originating from the main pulse in¬ 
teraction. 

The effect of the prepulse was modelled using the 
FLASH hydrodynamics solver in 2D cylindrical geometry. 
Helium was initialized following the experimental gas jet 
profile with uhc = 0.8 x 10^^ cm“^. The prepulse energy 
was modelled by depositing thermal energy evenly into 
a cylinder of gas at the laser focus (length 600/im and 
radius 80/im) centered near the initial critical density 
surface; riHe = Within Ins, a collisional blast wave 
is formed and expands out into the cold gas (HJ |23] . The 
blast wave forms a cavity with low density but high elec¬ 
tron temperatures at the center, surrounded by a cold 
but high density shock moving outwards. Fig.[^ shows 
the electron density at 25 ns after the prepulse, calculated 
from the Saha equation, for an initial 3 mJ of energy ab¬ 
sorbed in the plasma. Only the center of the cavity is 
hot enough to be ionized and detectable by interferom¬ 
etry, in agreement with the experimental measurements 
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FIG. 1. Electron plasma density from He targets, rie from 
interferometry: (a) immediately before, and (b) 300 ps after 
the LPI, for Epp ^ 220 mJ. (c) Ue and (d) total helium atomic 
density (irrespective of ionization), riHe from hydro simula¬ 
tions after 25 ns of expansion for 3 mJ absorbed energy. 


prior to the LPI in fig. Inspection of the neutral he¬ 
lium density (fig.[^) reveals the blast wave, which has a 
cavity wall of peak density given by that of a strong shock 
(7 ± l )/(7 — l)ni, where rii is the ion density and 7 , the 
ratio of specific heats (5/3 for He, 1.4 for H 2 )[ 2 T]. There¬ 
fore, using hydrogen with comparable laser parameters 
gives a higher density enhancement as well as a reduced 
scale length [22]. Even though this neutral profile exists 
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FIG. 2. Interferograms 300 ps after the LPI with He; (a) 
without prepulse, (b) with prepulse ideal for ion acceleration, 
Epp ^ 150 mJ, (c) and with prepulse too large for ion accel¬ 
eration, Epp ^ 1.27J. (d-f) corresponding nejncr- 
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FIG. 3. (a) Average flux of accelerated helium ions (/sr) 

> 0.1 MeV as sampled by the spectrometer (color scale) as 
a function of prepulse energy, Epp and main pulse ao. The 
white region indicates laser parameters tested experimentally, 
(b) Proton spectra for a H 2 target with Epp ^ 4 and 23 mJ, 
and He^ spectrum for a He target with Epp ^ 150 mJ (solid 
lines). The dotted lines are the detection thresholds, and the 
dashed red line is the deconvolved spectrum for Epp ^ 4mJ. 

before the LPI, it is only due to ionization by the fast 
electrons generated during the main interaction that it 
can be directly observed. Once ionized, the blast wave 
provides the ideal small scale-length plasma in which to 
generate collisionless shock waves. 

Interferograms for three different prepulse levels with 
helium gas are shown in fig.[^-c, with corresponding elec¬ 
tron densities in fig.[^-f. The images are all taken 300 ps 
after the LPI to reveal the full blast wave structure. Sim¬ 
ilar images are obtained with H 2 . To account for the ver¬ 
tical density gradient in the initial gas density, azimuthal 
symmetry is assumed in each half-cylinder above and be¬ 
low the laser axis. These regions were processed individ¬ 
ually by Abel inversion, with continuity assumed at their 
interface. 

With no prepulse (fig.[^ a & d), a small cavity of di¬ 
ameter ^ 100 /im forms around the laser focal position. 
The intense laser has self-focussed and channelled part 
way into the long density scale length (I 1 mm), cou¬ 
pling most of the energy into the plasma ramp. In this 
case, no forward accelerated ions were observed. Intro¬ 
ducing a prepulse with energy Epp ^ 150 mJ (fig.[^ b 
& e) shows a significant difference in the plasma distri¬ 
bution. A blast wave has been generated, so that the 
intense pulse interacts with a profile with steeper density 
gradient, ^100 /im, and higher peak density. Energetic 
ions were consistently generated in this regime. For fur¬ 
ther increase in Epp (fig.j^c & f), the prepulse induced 
blast wave propagates deeper into the jet and ultimately 
through it, leaving only an underdense remnant along 
the laser axis. For this higher prepulse, no ion beam was 
observed. Probing in this case, showed little change to 
the shape of the plasma before and after the LPI, except 
for a doubling of the density in the walls, consistent with 
ionization from He+ to He^+. 

Control of the helium density profile generated by the 
blast wave is essential for reproducible ion acceleration. 


Fig# shows the flux of detected He+ ions (color scale) 
as a function of main pulse normalized vector poten¬ 
tial, ao, and prepulse energy, Epp. Accelerated ions were 
only observed over a narrow range of prepulse energies, 
110 mJ < Epp < 220 mJ, showing the importance of op¬ 
timizing the density profile as in fig. [^b & d. Within this 
range of Epp and a^ ~ 1.4, the total flux of accelerated 
ions remains stable shot-to-shot at ^ 10 ^ ions/sr. 

Fig.j^ shows a He+ spectra obtained with Epp = 
150 mJ and a main pulse ao = 1.4, resulting in a broad 
energy distribution up to 1.5 MeV. This spectrum was 
typical for all the observed beams. Only He+ ions are 
observed although the main pulse laser intensity exceeds 
the threshold value to double ionize helium via field ion¬ 
ization (/ > 9x 10^^ Wcm“^). The observed single charge 
state is the result of charge transfer as the ions traverse 
the plasma and un-ionized gas [El ED, and recombina¬ 
tion as the accelerated ion beam co-moves with a low 
temperature electron cloud to the detector. 

For fixed rii and absorbed energy Eahsi a blast wave 
in hydrogen will give a steeper density ramp than for 
helium due to the lower 7 [ 22 ]. Spectra from shots 
with H 2 for different prepulse levels Epp = 4mJ and 
23 mJ are also shown in fig. Note that a signifi¬ 
cantly smaller prepulse was required for ion generation 
with hydrogen; observation of the relative blast wave 
size implies a significantly higher energy absorption. For 
Epp = 4mJ, a quasi-monoenergetic beam is observed 
with peak energy ^ 1.2 MeV and a deconvolved energy 
spread AEjE « 5 %, compared with a broadband beam 
with energies up to ~ 0.5 MeV for Epp = 23 mJ. As the 
radius for a spherical blast-wave oc {Eahs/'^iY^^ [H], 
lower Epp results in less expansion, providing a steeper 
target density gradient. Not only does this technique 
allow the production of reproducible shock accelerated 
ions, but for sufficiently steep gradients spectral shaping 
is also achievable. 

The LPI was modelled using the 2D PIC code EPOCH. 
The plasma was initialized with 30 particles per cell of 
cold ions (He^+ or H+) and electrons, in cells of size A/50. 
The laser had ao = 1.4, wq = 65/rm at focus, pulse 
length T = 5ps (FWHM), and was linearly polarized 
in the transverse plane. The initial density distribution 
was varied to investigate the effect of the optical density 
tailoring. 

To model the LPI with no prepulse, the initial density 
profile was set to rise linearly from 0 to 2 Uc over 800 /im 
with He^+. The laser self-focuses and bores a channel 
in the plasma, shown in fig. All the laser energy is 
expended in the underdense region and no forward shock 
is formed, although lower velocity transverse electrostatic 
shocks are observed |3]. 

The effect of the prepulse was simulated using a plasma 
density distribution extracted from the interferometry for 
Epp ^ 150 mJ (green box in fig. #). The profile was 
taken along a line radially off-set from the laser axis in 
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order to reduce errors introduced by the Abel inversion. 
In the blast-wave ramp the density is well approximated 
by n = where Umax ~ 6ric is the peak 

density, zq is the position of the peak and lexp ~ 40 /im is 
the scale length. In this simulation, the laser pulse shows 
a modest increase in intensity due to self-focussing in the 
ramp, but the short ramp length means the laser is not 
adversely absorbed. The pulse penetrates up to the criti¬ 
cal surface where a combination of radiation and thermal 
pressure launches a collisionless shock (fig. &) . Ions are 
reflected at the shock front to twice the shock speed 
The result agrees well with optical shadowgraphy 300 ps 
after the LPI (inset of fig.[^), which shows remnants 
of the accelerating collisionless shock emerging from the 
collisional prepulse induced blast wave. 

The z — Pz phase space (fig.[^) demonstrates that the 
ions originate from the position of the shock, but also 
that the ions are reflected with a large energy spread. 
Fig.|^ shows the position of the shock front and critical 
surface on the laser-axis as a function of time for this sim¬ 
ulation. The initial shock velocity Vgh ~ 5 x 10^ ms“^ ex¬ 
ceeds the measured hole-boring velocity Vhb = 3.1 x 10^ ~ 

^LI Haim- Ions reflected off the shock at its 

peak velocity would gain ^3.1 MeV, in agreement with 
simulated maximum ion energies of ^ 3.3 MeV (fig.[^). 
Scaling the simulated ion spectra to the experimental 
spectra and imposing the experimental detection limit 
as in fig. would result in a maximum detectable ion 
energy of ^ 2 MeV, consistent with that measured. The 
energies in the case of a shaped target are significantly 
larger than for the no prepulse case. 

The shock slows at the end of the laser drive while still 
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FIG. 4. Results of PIC simulations. Ion density map 12 ps 
after LPI for: (a) no prepulse, (b) experimental initial profile 
for Epp = 150 mJ obtained from fig. §• (inset is optical shad¬ 
owgraphy showing shock remnants) and (c) z-pz phase space, 
(d) The z-pz phase space for for hydrogen plasma with expo¬ 
nential ramp {lexp = 20 pm.) 
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FIG. 5. (a) Critical surface (red) and shock (blue) position 

with time for simulation shown in figJ^Vc, and shock po¬ 
sition for hydrogen simulation in fig.ppT (b) Comparison of 
He^^ spectra for 800 pm linear ramp (no prepulse) (green), 
prepulse fig.[^Vb (blue), and for hydrogen with exponential 
profile, lexp = 20pm (black). 


reflecting ions, resulting in the large momentum spread in 
fig.|4)3 and the corresponding broadband axial ion spec¬ 
trum in fig. HI- The shock deceleration is mainly due 
to spherical expansion of the front. However, the ex¬ 
periment showed spectrally peaked proton beams were 
produced from sharper density gradients using hydrogen. 
Density profiles with an exponential ramp, and variable 
length scale, were used to reproduce the sharper profiles 
achievable with hydrogen. At the edge of the blast-wave, 
the density dropped to Umax I the density drop for a 
strong shock, imitating the experimental profile. 

A scale length lexp = 50 pm produces a broadband 
beam, as described for helium. However, shortening the 
scale length to lexp = 20 pm results in the shock breaking 
through the density discontinuity at the blast-wave front. 
The sudden change in thermal pressure ahead of the front 
triggers the shock immediate dissipation, just at the end 
of the LPI (fig.[^). The axial z —pz phase space, fig.[^, 
demonstrates the formation of a single bunch with veloc¬ 
ity near 2vs ~ 13 x 10^ m/s. Particle reflection stops and 
the peaked spectrum is maintained (fig.|^). This gen¬ 
eration of peaked spectra differs from thermally driven 
shocks in isothermal plasmas 0E], where reflection at a 
uniform speed persists over a longer time scale. Further 
simulation with lexp < 10 pm with helium also demon¬ 
strated spectrally peaked beams. Producing such gra¬ 
dients in helium would be experimentally achievable by 
shortening the time between the prepulse and the drive 
pulse, which was not possible with our experimental ge¬ 
ometry. 

In this work, we have observed shock acceleration of 
protons and helium ions. The initial plasma density pro¬ 
file is found to be critical to achieve not only shock accel¬ 
eration but also spectral control of the ion beam. Using 
hydrogen, we found that long density gradients lead to 
the production of broadband shock accelerated beams, 
while steeper gradients allow for the generation of quasi 
monoenergetic beams. Helium ions have also been ac¬ 
celerated using this scheme, demonstrating that it could 
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be extended to other high Z gaseous species, providing a 
route towards a wide variety of easily replenishable high- 
repetition rate ion sources for future nuclear physics ap¬ 
plications. 
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